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ABSTRACT 


The  purpose  of  this  report  is  to  present  a  breakdown  surface  model  for 
evaluating  thermal  backscattering  flux  from  the  supersonic  exhaust  plume  of  a 
space-based  HF  laser.  The  plume  is  of  ring  symmetry.  It  consists  of  a 
gaseous  mixture  of  H,  HF,  H2 ,  DF  and  He.  Fluxes  of  these  species  are 
considered  separately.  The  model  is  carefully  analyzed  and  is  shown  to 
overestimate  the  flux.  Actual  flux  levels  of  the  heavy  corrosive  molecules 
(HF,  DF)  have  been  found  to  be  exceedingly  low.  It  is  concluded  that  the 
contribution  of  thermal  backscattering  to  contaminating  flux  of  HF  and  DF  can 
be  neglected.  This  work  is  an  extension  and  modification  of  the  recent  thesis 
work  done  by  S.  E,  McCarty  at  the  Naval  Postgraduate  School. 
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1.  INTRODUCTION 


This  report  is  a  presentation  of  one  part  of  a  study  on  the  contaminating 
backflow  from  the  exhaust  plume  of  a  large  space-based  HF  laser  (Figure  1). 
The  exhaust  plume  is  an  underexpanded  supersonic  ring-jet,  designed  to  stay 
clear  of  the  spacecraft  by  maintaining  a  Prandt 1-Meyer  turning  angle  at  the 
nozzle  lips  of  well  below  90°.  However,  it  is  wpll  known  from  experience  with 
rocket  plumes  in  space  11*2]  that  cavity  regions  (where  continuum  gasdynamic 
theory  predicts  vacuum)  are  filled  with  a  free-molecular  flow.  This  back  flow 
is  largely  due  to  viscous  effects,  which  give  rise  to  a  "spill-over"  of  the 
boundary  layer  around  the  nozzle  lip (3].  Assuming  the  boundary  layer  can  be 
eliminated  (e.g.,  an  expanding  step  design  of  the  nozzle  lip),  there  are 
two  more  mechanisms  which  lead  to  backflow:  thermal  backscatter ing  and 

scattering  by  ambient  molecules  traveling  at  orbital  speeds.  Since  these 
effects  are  a  small  perturbation  to  the  exhaust  flow  field,  they  can  be 
considered  independently  (the  total  backflow  will  be  a  superposition  of 
contributions  due  to  these  two  effects).  As  a  first  phase  of  our  broader 
study,  we  consider  solely  the  contribution  of  thermal  backscattering  to  back- 
flow  from  a  ring-plume  of  an  HF  laser,  via  a  simple  model  of  molecular 
effusion  from  a  breakdown  surface,  fashioned  after  ideas  suggested  by 
NollerlM.  Our  results  indicate  that  the  backflow  of  the  heavier  contaminants 
(HF,  DF)  due  to  thermal  backscattering  is  negligible. 

Naturally,  our  study  pertains  to  presumably  typical  operating  conditions 
of  the  HF  laser.  These  operating  conditions  were  largely  determined  from  a 
report  on  some  HF  laser  tests  conducted  at  TRW  in  1971 [5]  (in  particular. 
Table  5,  Test  III,  of  this  report).  The  typical  parameters  at  the  nozzle  exit 
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are: 

Composition  (mole  fractions): 

Specific  Heats  Ratio: 

Mach  Numbe 

Average  Molecular  Weight 
Stagnation  Temperature 
Stagnation  Density 
Molecular  Diameter,  assuming  it  is 
uniform  for  all  species  (hard- 
sphere  collisions) 


[HJ  -  .091,  [HF]  -  .091, 

[H21  -  .104,  [DF]  =  .135, 

[He)  -  .579 

Y  *  1.54  (assuming  ideal  gas) 
Mi  ■  4.0 

=  7.27  [kg/kg  mole]  (1.1) 
T0  -  2300  [K] 
p0  -  0.0075  [kg/m3] 

D  *  2.5  x  10"10  [m] 


The  exit  Mach  number  can  be  chosen  higher  than  Mi  ■  4,  but  not 
considerably  lower  than  this  value,  since  Mi  *  4  results  in  a  modest  clearance 
angle  of  41°  between  the  limiting  (vacuum)  characteristic  of  the  lip-centered 
rarefaction  fan  and  the  spacecraft.  We  assume  isentropic  flow  throughout  the 
diffuser so  that  upon  specifying  the  composition  and  flow  variables  at  the 
diffuser  inlet,  along  with  Mi  at  the  diffuser  exit,  the  exit  flow  is  fully 
determined.  One  exception  to  this  definition,  however,  is  the  stagnation 
temperature,  which  was  estimated  as  T0  =  1400  [K]  at  the  diffuser  inlet t ^ ] . 
We  set  T0  ■  2300  [K],  which  corresponds  to  complete  hydrogen  recombination, 
even  though  the  flow  in  the  diffuser  is  of  a  nearly  frozen  composition  due  to 
the  low  rate  of  hydrogen  recombination t 3 ] .  The  reason  for  this  choice  is  that 
given  the  uncertainty  in  determining  T0,  which  results  from  an  uncertainty  in 
the  degree  of  hydrogen  recombination,  it  is  the  most  pessimistic  choice, 
resulting  in  higher  thermally  backscattered  flux. 

The  model  that  we  propose  for  evaluating  the  backscattered  flux  arriving 
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at  the  spacecraft  (Figures  1,  2)  is  based  on  the  effusive  breakdown  surface 
concept  suggested  by  Nollerf^.  The  gradual  transition  from  continuum  to 
collisionless  flow,  which  invariably  takes  place  at  the  outer  fringes  of 
exhaust  plumes  having  a  near-vacuum  background  environment,  is  replaced  by  an 
abrupt  change.  We  assume  that  the  flow  regime  in  each  stream  tube  changes 
from  continuum  (with  local  thermodynnamic  equilibrium)  to  collisionless,  upon 
crossing  some  breakdown  surface. 

An  important  simplification  is  introduced  in  the  case  of  a  large-radius 
spacecraft  (about  2 . 5 [ m ] ) ,  by  observing  that  the  temperature  along  the 
breakdown  surface  decreases  so-  sharply  with  the  distance  from  the  nozzle  lip, 
that  the  segment  contributing  significantly  to  thermal  backscattering  is  only 
about  0.01  to  0.1 [m]  long.  Consequently,  the  lip-centered  rarefaction 
ring-fan  may  well  be  approximated  by  the  standard  (planar)  Prandtl-Meyer 
flow  field. 

The  structure  of  this  report  is  as  follows.  The  breakdown  surface  and  the 
molecular  effusion  flux  from  it  are  obtained  in  closed-form  expressions  in 
Section  2.  Section  3  is  devoted  to  the  spatial  integration  scheme,  which  is 
the  evaluation  of  the  flux  arriving  at  a  certain  point  on  the  spacecraft. 
Results  of  flux  distribution  along  the  spacecraft  for  the  presumed  laser 
operating  range  are  presented  and  discussed  in  Section  4,  followed  by  a 
critical  examination  of  the  breakdown  surface  model.  Conclusions  are  given 
in  Section  5,  and  Section  6  is  a  list  of  references.  The  code  RINGBD,  which 
computes  the  flux  by  numerical  integration  over  the  breakdown  (effusing) 
surface,  is  given  in  Appendix  A. 
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2.  BREAKDOWN  SURFACE  AND  EFFUSION  FLUX 


Our  model  for  the  thermally  backscattered  flux  arriving  at  the  surface  of 
the  spacecraft  is  essentially  a  modification  of  Noller's  concept  of  a 
breakdown  effusive  surface We  substitute  his  definition  of  a  breakdown 
surface  by  a  similar  one  introduced  by  Birdf^*  Section  8.3],  We  obtain  the 
one-sided  effusion  flux  from  the  breakdown  surface  by  integrating  over 
velocity  space  as  suggested  by  Noller^^,  except  for  the  fact  that  we  compute 
flux  rather  than  density  and  we  also  consider  the  flux  of  species  having 
molecular  weight  different  from  the  average.  In  the  following,  each  one  of 
these  steps  is  described  in  some  detail,  beginning  with  the  breakdown  surface. 

As  mentioned  in  the  introduction,  the  lip-centered  rarefaction  fan  is 
approximated  by  a  planar  Prandtl-Meyer  flow  field  (Figure  2).  The  standard 
expressions  for  this  flow  field  have  Mach  number  (M)  as  the  independent 
parameter,  thus  M  varies  between  M  =  Mj  at  the  exit  and  H  +  «  at  the  limiting 
(vacuum)  characteristic.  (Index  1  always  refers  to  exit  conditions,  i.e.,  to 
parameters  evaluated  at  M  =  M^. 

*(H)  =  -  ?(M)  +  ?1  +  vf  f 

C(M)  =  YiU  ARCTAN[r-l/2(M2-l)W2j  r=(-iiy) 

U(H]  =  ARCS IN( M-1 )  (2.1) 

e(M)  -  <|(M)  -  p(M) 

where  y  is  the  angle  of  characteristic  lines,  and  9  is  the  angle  of  the 
velocity  vector  (or  streamline). 
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Adopting  Bird's  definition  of  a  breakdown  parameter,  which  was  first 
introduced  in  conjunction  with  a  spherical  source  flow^  »  Section  8.3]  ancj 
later  was  shown  to  be  meaningful  also  in  a  Prand tl-Meyer  flow!'7},  we  define 
the  breakdown  surface  as  having  a  constant  value  of  B,  where  B  is  given  by: 


B 


y  i  dp 
v  p  dS 


(2.2) 


Here  p,  U,  v,  S  are  local  flow  density,  speed,  collision  frequency,  and 
coordinate  along  streamlines  (thus  restricting  this  definition  of  B  to 
stationary  flows).  Prom  the  geometrical  relationships  in  a  Prand tl-Meyer  fan 
(Figure  2)  and  from  (2.1)  we  get: 

7s  =  “  (1/R)  ^  sin  w  =  -  ‘  1)1/2^  (p/R) 

1  (2.3) 

P(M)  =  Po  (1  +  M2)  *-l 

Using  the  expression  for  collision  frequency!^!; 


>0  =  ^n)uz  (nod2c0; 


(2.4) 


where  N0,CQ,D  are  stagnation  number  density,  stagnation  sound  speed,  molecular 
diameter,  and  using  U  =  MC  in  conjunction  with  (2.2)  and  (2.3),  we  get: 


Rb(M)  =  (BNqD  ) 


2i-i  mi 


1/2 


1 

;m2-i)1/2  (i+  m2)  t-1 


(2.5) 


2l^rr  1 

This  expression  is  almost  identical  to  that  of  Bi rd ! ^ I ,  the  main  difference 


being  in  assuming  a  constant  collision  diameter  (hard  spheres),  which  we 


believe  to  be  commensurate  with  the  overall  crudeness  of  the  model.  The 
breakdown  surface  as  defined  by  equations  2.5  and  2.1,  starts  at  point 
[RgCH^),  on  the  exit  characteristic  M  ■  M  ,  which  we  refer  to  as  the 

intial  point  (see  Figure  4).  However,  a  breakdown  in  continuum  flow  also 
takes  place  on  the  segment  of  the  exit  characteristic  between  the  corner  and 
the  intial  point,  since  the  value  of  the  breakdown  parameter  there  (Equation 
2.2)  is  clearly  larger  than  the  value  of  B  ufced  in  defining  the  breakdown 
surface  (Equation  2.5).  Hence,  the  breakdown  surface  defined  by  (2.5)  should 
be  supplemented  by  that  segment.  We  refer  to  the  combined  surface  as  the 
augmented  breakdown  surface.  The  segment  on  the  exit  characteristic  is 
referred  to  as  the  supplementary  breakdown  surface. 


The  one-sided  directed  effusion  flux  is  defined  as  the  number  flux  of 
molecules  per  unit  area  of  an  area  element  normal  to  the  flux  direction,  per 
unit  solid  angle  about  the  flux  direction.  It  is  obtained  as  a  function  of 
local  Mach  number  and  the  angle  <  between  the  flux  direction  and  the  local 
velocity  vector,  by  repeating  Noller's  velocity  integration  scheme^’^*  (6)1, 
with  an  added  factor  of  molecular  speed  in  order  to  obtain  flux  (rather  than 
density  as  in  Noller's  work).  The  resulting  expression  for  species  i  is 
readily  obtained  by  using  standard  definite  integrals: 


Fi(fi)  -  h*  N0CQ  (WA/W±)l/2  (1  +1ZL  M2)  2<^_1) 

[(2yit3)~1/2  (l  +  (1/2)y!12cos2<)  EXP  (-  yM2/2)  + 

( 2 tt )"~ 1  M  (3/2  +  (1/2)yM2cos2<)  cos<  EXP  (-  (1/2)yM2  sin2*) 


(2.6) 


ERFC  (-  (Y 12)  M  cos*)] 


CO 

ERFC(v)  =  2  n-1/2  f  EXP  (-x2)dx  (complementary  error  function) 

v 
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(2.6  Continued) 


M  -  (Wt/WA)1/2  M 

-  Mole  fraction  of  species  i. 

Wi  -  Molecular  weight  of  species  i. 

The  dependence  of  F^(M)  on  the  flux  angle  k  is  so  sensitive  that  for  some 
Mach  number  around  M  *  10,  the  backflow  (in  the  typical  operating  range)  is 
virtually  negligible.  In  the  following  section  we  describe  how  the  flux  F^(M) 
is  integrated  over  the  augmented  breakdown  surface,  yielding  the  backscattered 
flux  arriving  at  the  surface  of  the  spacecraft. 


3.  FLUX  INTEGRATION 


The  effusion  flux  Fi(M)  given  by  (2.6)  above,  is  defined  in  such  a  way 
that  the  number  of  molecules  effusing  from  an  area  element  AAg  of  the  break¬ 
down  surface  and  arriving  at  an  area  element  AAg  on  the  spacecraft  (per 
second),  is  given  by: 

_2 

Fi(tl)  (AAg  COS  ag)  (AAg  COSag)  L  [molecules  per  second]  (3.1) 

u  o 

where  ag,  ag  are  the  angles  between  the  line-of -sight  Lgg  (Figure  3)  and  the 
normals  to  the  breakdown  surface  and  the  spacecraft  surface  respectively.  Lgg 
is  the  distance  between  the  elements  AAg  and  AAg.  Dividing  equation  (3.1)  by 
AAS  and  integrating  over  the  breakdown  surface,  the  flux  per  unit  area  of 
the  spacecraft  is  given  by: 

Qi  ■  jF^(M)COSClg  COSOg  L~j;  dAg  [Molecules  per  second  per  m^j  (3.2) 

D  j 

The  integration  scheme  for  Qj  over  the  breakdown  surface  is  expressed  in 
terms  of  the  set  of  polar  coordinates  R,  \|*,  $  (Figure  3).  For  a  point 
(R,  i p,  4>)  on  the  breakdown  surface,  using  Cartesian  coordinates  (X,  Y,  Z)  and 
the  angle  u  between  X-axis  and  the  line-of-sight  Lgg  we  obtain  the  following 
geometrical  relationships: 

X  -  R  COS«  ;  Y  ■  (A0  +  R  Silty)  COS$  ;  Z  «=  (A0  +  R  SIN<p)SIN<p 

C0SWM)  "  (ao+R<M)SIN\jJ  (3.3) 

U  -  U(COS0,  SIN0  COS«p ,  SIN0  SIN$) 

Lbs/I^BsI  “  (COSui,  SINo)  COSB,  SINw  SINB) 

Z 

(Y-A0) 


TAN  8 
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The  cosines  COS<,  COSag,  COSag  in  (3.2;  are  expressed  as  scalar  products 
of  [Lbs/ I LBS 1 1  ancl  unit  vectors  along  the  local  velocity  vector  U,  along 
the  local  normal  to  the  breakdown  surface  and  along  the  local  normal  to  the 
spacecraft  surface,  correspondingly. 

The  integration  is  performed  numerically  in  two  phases,  the  first  being 
the  integration  along  the  supplimentary  breakdown  surface  (Figure  4).  For 
this  first  phase,  the  straight  line  segment  which  constitutes  the 
supplementary  breakdown  surface  is  divided  into  several  intervals  of  length  AR 
(typically  10  intervals).  Each  interval  generates  a  half-strip  by  rotating  it 
from  $  *  0  through  <j>  =  )»  This  strip  is  in  turn  subdivided  into 
several  sub-intervals  of  A<J>  each  (typically  10  intervals).  The  total  flux 
arriving  at  Xs  is  obtained  by  summing  contributions  from  each  sub-interval 
(two-dimensional  integration).  When  the  integration  along  the  supplementary 
breakdown  surface  is  concluded,  it  is  continued  into  the  breakdown  surface, 
where  AR  intervals  are  replaced  by  breakdown  surface  intervals  that  correspond 
to  a  fixed  Mach  number  increment  AM  (typically  AM=0.1).  The  integration 
proceeds  along  the  breakdown  surface  (2.5)  until  the  contribution  of  the  last 
AM  strip  is  negligibly  small.  The  computation  time  is  modest  (about  1  second 
CPU  per  Xs  point,  on  IBM  3033  mainframe  computer).  The  computations  were 
carried  out  by  a  code  RINGBD  written  specifically  for  this  purpose.  Further 
details  of  the  scheme  and  programming  can  be  obtained  by  reading  this  code 
which  is  given  in  Appendix  A. 
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4.  RESULTS  AND  DISCUSSION 

4.1  Presentation  of  Results 

The  molecular  flux  backscattered  to  the  spacecraft  from  the  surface 
of  continuum  flow  breakdown  in  the  lip-centered  rarefaction  fan,  has  been 
computed  for  all  five  species  H,  HF,  H2,  DF,  He.  The  results  are  depicted  in 
Figures  5  to  9  respectively.  For  each  species  two  more  cases  were  computed  in 
addition  to  the  nominal  case  (1.1),  where  the  stagnation  density  p0  was 
replaced  by  po/10  and  by  po*10  (see  Figures  5  to  9).  This  has  been  done  in 
order  to  demonstrate  the  effect  of  variations  in  exit  flow  conditions  on  the 
flux.  The  particular  choice  of  p0  was  motivated  by  the  fact  that  the  effects 
of  changing  pD  are  not  obvious.  The  effects  of  changing  the  exit  Mach  number 
or  the  stagnation  temperature  T0  are  rather  obvious  (a  higher  flux  would 
result  from  either  a  decrease  in  11  or  an  increase  in  T0).  It  turns  out  that 

for  points  lying  not  too  near  the  nozzle  lip  (Xs  >  0.1  m) ,  the  lower  density 

flow  generates  a  higher  backscattered  flux! 

In  addition  to  varying  p0,  we  also  varied  the  breakdown  parameter  B, 
obtaining  a  surprising  result.  The  computation  was  performed  for  a  particular 
species  (HF),  and  the  results  obtained  upon  replacing  B=0.05  (nominal  value) 
by  B/2  and  by  B*2  are  brought  in  Figure  10. 

It  turns  out  that  the  B/2  case  has  the  higher  flux.  This  is  somewhat 
surprising,  since  a  lower  value  of  B  in  a  centered  rarefaction  fan  (equation 
2.5)  means  that  the  breakdown  of  continuum  flow  takes  place  in  a  region 
further  out  from  the  corner.  In  a  source  flow  (e.g.,  a  spherical  source), 

that  implies  lower  density  and  temperature,  which  would  give  rise  to 


lower  thermally  backscattered  flux 


An  explanation  to  these  seemingly  counterintuitive  results,  along 
with  some  deeper  insight  into  the  breakdown  surface  model  as  it  is  applied  to 
a  centered  rarefaction  flow,  can  be  obtained  by  taking  a  close  look  at  the 
flow  field  and  the  breakdown  surface  in  the  vicinity  of  the  corner.  We  take 
up  this  matter  in  the  following  sections. 

We  conclude  the  presentation  of  results,  by  comparing  the  flux  (in 
the  nominal  case)  of  the  five  species  with  each  other  (Figure  11).  This 
figure  underlines  the  fact  that  the  flux  of  light  species  (H,  H2 »  He)  is  many 
orders  of  magnitude  (typically  10 15)  times  that  of  heavy  species  (HF,  DF). 
Indeed,  these  results  demonstrate  a  well  known  effect:  When  an  expanding 

gaseous  mixture  of  light  and  heavy  molecules  experiences  a  breakdown  of 
continuum  flow,  a  separation  of  species  takes  place  (see  e.g.,  the  work  of 
Cattolica  et.  al.  [8]). 

4.2  The  Breakdown  Surface  and  Streamlines 

Consider  the  parametric  description  Rg(M)  for  the  breakdown  surface 
(Equation  2.5).  Normalizing  R  relative  to  the  exit  mean  free  path  Xj,  we  get: 

1 

Rg(H)  -  RB(M1)[(M2  -  1)/(M*  -  1)]1/2  [(1  +^li  M*j/(l  il2)]Y_1 

M[1i)Ai  *  [(W2)1/2/(rt-i)B]  (m2  -  ij1/2 

(4.1) 

1 

,  1/2  “ 1  r  v_i  2  v„.i 

*1  “  l2  *D2N0)  [1  +Xl  Mj  * 

The  normalized  surface  Rg(M)/X^  is  thus  independent  of  stagnation  density, 
depending  only  on  y,  M  ,  B. 
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Let  us  now  derive  a  parametric  equation  Rg(M)  for  a  streamline  that 
enters  the  fan  at  point  RgCM^  on  the  exit  characteristic.  The  following 
geometrical  relationship  is  readily  obtained  by  considering  two  characteristic 


lines  ^  and  ip-t-Aip  and  a  streamline  inclined  at  the  Mach  angle  p  to  them: 


dRs(<|>) 


-  RgCty)  (tanp)'1 


(4.2) 


Using  the  standard  Prandtl-Meyer  functions  (2.1),  we  get  the 


following  differential  equation  for  Rg(M): 


_ *  -ISsi rc ).  (Xtl)  i  +  m2  )_1 

r„(M)  dM  1  2  '  "l1  +  2  "  j 


(4.3) 


This  equation  is  readily  integrated,  giving: 

_I±i _ 

RS(M)  -  R^M,)  [(1  +  IliM2)/(l  +  2(Y‘°  (4.4) 

As  pointed  out  by  tiird^l,  there  is  a  particular  streamline  Rga(M) 
which  asymptotically  approaches  the  breakdown  surface  for  large  M,  since  the 
ratio  R8(M)/RB(M)  tends  to  a  constant  (not  zero)  when  M  *  ®,  (Strictly 
speaking,  this  holds  only  for  hard-sphere  molecules,  i.e.,  only  when  u>  -  0.3 
in t ^ 1 ) .  The  limit  is: 


lim  MM)  Rs(Hi>  f  Mf  -1  il/2 
Rb(m5  ”  ^('h)  ^2  +  2  j 

1  y-1 


(4.5) 


For  the  limiting  streamline  Rga(M) ,  the  ratio  Rsa(M)/Rg(M)  should 
tend  to  1.  This  determines  the  point  Rggdij)  at  which  the  limiting  streamline 
enters  the  fan,  as  well  as  the  entire  line  Kga(M): 


2  2 

Ml  +  -=7,1/2 


R8a(M1)  -  RB(Mt)  [  ---gT-j  ] 

M.  -  1 


(4.6) 
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Clearly,  Rsa(M)  is  larger  than  Rg(M)  for  any  M  >  Mj,  so  that  no 
streamline  beyond  Rsa(M)  can  cross  the  breakdown  surface.  This  pattern  is 
shown  in  Figure  4,  where  Rga(M)  is  denoted  "streamline  2",  and  "streamline  1" 
is  the  streamline  R^M,)  =  RbCMj). 

All  this  leads  to  the  following  observation  regarding  the  continuum 
breakdown  of  the  flow  in  a  centered  rarefaction  fan^l.  Referring  to  Figure 
4,  the  fluid  entering  the  fan  through  the  supplementary  breakdown  surface 
(i.e.,  through  the  exit  characteristic  between  the  corner  and  the  initial 
point  of  streamline  1),  experiences  breakdown  immediately  upon  crossing  this 
surface.  Every  streamline  between  streamline  1  and  streamline  2  crosses  the 
breakdown  surface  at  some  Mach  number  M  >  ,  and  at  that  point  the 
continuum  flow  breaks  down.  All  fluid  entering  the  fan  beyond  stream-  line  2 
will  never  pass  through  the  breakdown  surface,  and  hence  will  maintain  a 
continuum  flow  regime  all  the  way  to  infinity.  Of  course,  that  is  only  true 
for  planar  centered  rarefaction  fans.  When  the  exhaust  flow  emerges  from  a 
nozzle  of  finite  width,  and  especially  when  the  exhaust  jet  has  a  ring 
symmetry  (as  in  our  case),  the  breakdown  surface  gradually  curves  in  a 
balloonlike  shape  towards  the  opposite  nozzle  lip,  forming  the  familiar  plume 
pattern  (Figure  1). 

4 . 3  Analysis  and  Discussion  of  Results 

The  foregoing  analysis  is  now  used  to  explain  the  variation  in  back- 
scattered  flux  due  to  a  change  in  exhaust  flow  conditions  at  the  nozzle  exit. 
Specifically,  we  consider  a  tenfold  decrease  in  stagnation  density  (i.e.,  the 
case  po/10),  and  hence  a  tenfold  increase  in  the  exit  mean  free  path  X  , 

The  effusion  flux  from  the  breakdown  surface  is  proportional  to  the 
local  density,  so  one  would  expect  to  observe  a  decrease  in  flux,  rather  than 


3 


an  increase  (see  Figures  5  to  9,  for  Xs  >0.1  m) .  Other  factors  causing 
increased  flux,  must  then  be  larger  than  10  so  that  they  more  than  offset  the 
1/10  factor  in  density.  It  turns  out  that  these  effects  are  mainly 
geometrical,  in  that  a  tenfold  increase  in  X^  causes  the  domain  of  integration 
on  the  breakdown  surface  to  increase  more  than  tenfold.  In  the  (X,Y)  plane 
there  is  a  tenfold  "blowup"  of  the  breakdown  surface,  due  to  the  self-similar 
structure  of  the  Prandtl-Meyer  flow  field.  As  a  result  of  this  "blowup”  in 
(X , Y) ,  the  angular  integration  range  $max  also  increases,  albeit  not  linearly 
(Equation  3.3).  Another  geometrical  effect  is  an  increase  in  the  flux 
incidence  cosine  factor  cosas  (see  Equation  3.2),  which  for  points  Xg 
sufficiently  far  from  the  nozzle  lip,  increases  roughly  tenfold  (while  the 
other  cosine  factor  cosag  is  almost  constant).  All  this  provides  a 
qualitative  explanation  for  the  observed  increase  in  flux  at  far  points  (Xs  > 
0.1  m) . 


As  for  the  near  range  (Xs  <  0.1  m) ,  another  effect  becomes 
increasingly  significant  as  Xs  approaches  the  nozzle  lip.  The  turning  angle 
k,  by  which  backscattered  molecules  have  to  be  deflected  relative  to  the  flow 
velocity  vector  in  order  to  reach  point  Xs  on  the  spacecraft  (Figure  2), 
increases  with  the  size  of  the  breakdown  surface  (fixed  M  and  Xs).  Since  the 
local  effusion  flux  (Equation  2.6)  decreases  rather  sharply  as  <  is  increased, 
the  net  result  is  a  tendency  to  get  a  reduced  backscattered  flux  at  near 
points  such  as  Xg  “  0.01  m  (Figures  5  to  9). 

We  now  turn  to  the  effect  of  changing  the  value  of  the  breakdown 
parameter  B.  From  equation  4.1  it  is  clear  that  multiplying  B  by  some  factor 
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will  have  the  same  "blowup"  effect  as  dividing  M  by  the  same  factor.  A 

tenfold  decrease  in  B  is  thus  geometrically  equivalent  to  a  tenfold  decrease 

in  pQ.  However,  since  the  local  effusion  flux  at  the  breakdown  surface  is 
proportional  to  p0  while  it  is  independent  of  B,  the  B/10  case  will  have  ten 
times  as  much  backscattered  flux  as  the  pQ/ 10  case.  In  order  to  illustrate 
the  sensitivity  of  the  flux  estimates  to  an  uncertainty  in  the  appropriate 
value  of  B,  we  computed  the  cases  B/2  and  B*2  for  one  species  (HF),  and  they 
are  presented  in  Figure  10.  The  variation  in  flux  relative  to  the  nominal 
case  (B  =  0.05),  is  by  a  factor  no  larger  than  about  5.  Results  for  other 

species  were  found  to  exhibit  comparable  variations. 

Does  this  observation  about  the  dependence  of  the  breakdown  surface 

on  B  agree  with  the  breakdown  surface  appropriate  to  the  far  field  of  the 

exhaust  plume?  In  stationary  source  flow  into  vacuum,  and  when  M  »  1,  the 

6-1 

breakdown  parameter  varies  with  radius  as  B  ~  R  (5s*  1  for  cylindrical 
source,  6=2  for  spherical  source).  In  a  ringjet,  the  stream  tubes  of  the 
exhaust  plume  generally  diverge  at  a  rate  higher  than  that  of  stream  tubes  in 
a  cylindrical  source  flow,  so  the  effective  value  of  6  in  a  ringjet  is  5  >  1. 
Hence,  in  this  case  the  far  field  breakdown  surface  moves  downstream  along 
each  stream  tube  as  the  value  of  B  increases.  This  is  indeed  geometrically 

compatible  with  the  fact  that  near  the  corner  of  the  lip-centered  rarefaction 

_1 

fan  B  ~  R  ,  as  shown  schematically  in  Figure  12.  The  dependence  of  the  break¬ 
down  surface  on  B  near  the  corner  and  in  the  far  field,  thus  assures  that 
complete  breakdown  surfaces  corresponding  to  different  values  of  B,  do  not 
intersect  (Figure  12). 

In  the  foregoing  discussion  it  was  pointed  out  that  variations  in 
flux  caused  by  changes  in  parameters  such  as  p0  and  B,  were  directly  related 
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to  the  self-similar  structure  of  the  Prandtl-Meyer  flow  field.  It  has  been 
further  shown  that  these  variations  are  well-understood  within  the  framework 
of  the  breakdown  surface  model  and  that  they  are  not  excessively  large.  Are 
we  to  conclude  that  the  thermally  backscattered  flux  estimates  of  the  present 
model  are  also  physically  plausible  and  reliable?  In  the  following  section  we 
take  up  this  matter,  arriving  at  some  interesting  conclusions  about  this  model 


4.4  Cr i tlcal  Examination  of  the  Mode  1 

Consider  the  centered  rarefaction  flow  field  of  a  compressible  fluid 
negotiating  an  expansive  corner  at  supersonic  speed  (Prandtl-Meyer  flow).  The 
streamlines  of  this  flow  field  have  an  orderly  "layered"  structure,  with  each 
streamline  curving  around  the  corner,  starting  at  its  point  of  entrance  into 
the  fan  (see  Figure  4). 

The  present  model  is  based  on  the  stipulation  that  there  is  a  point 
of  continuum  flow  breakdown  on  each  streamline,  provided  this  streamline  is 
not  beyond  a  certain  limiting  streamline.  Consider  a  sample  molecule  effusing 
from  this  breakdown  point  toward  the  spacecraft.  It  advances  at  constant 
speed  along  a  straight  line  trajectory,  traversing  all  inner  streamlines. 
Since  the  flow  velocity  vector  points  away  from  the  spacecraft,  and  since  the 
flow  is  highly  supersonic  so  that  the  velocity  of  most  individual  molecules 
does  not  differ  much  from  the  flow  velocity  (i.e.,  it  is  a  "cold”  flow),  any 
collision  of  the  sample  molecule  with  a  mainflow  molecule  will  most  probably 
divert  the  sample  molecule  away  from  the  spacecraft..  What  is  the  probability 
that  a  sample  molecule  would  traverse  this  cross  flow  collisionlessly ?  This 
probability  is  simply  exp(-n),  where  n  is  the  expected  number  of  collisions 
along  the  straight-line  trajectory  from  the  point  of  breakdown  to  the 
spacecraft.  In  the  typical  operating  conditions  assumed  here,  we  estimated  n 
to  be  roughly  about  10.  Since  this  no-collision  probability  factor  is  ignored 
in  the  formulation  of  the  present  model,  the  backscattered  flux  may  be 
exaggerated  by  a  factor  of  exp(10)  or  about  104.  We  conclude  that  in  all 
likelihood,  the  prediction  of  the  breakdown  surface  model  for  thermally 
backscattered  flux  from  a  centered  rarefaction  flow,  is  substantially 


overestimated . 


Can  anything  be  done  to  improve  the  present  model?  One  may  be 
inclined  to  suggest  at  this  point  that  the  obvious  remedy  is  to  incorporate 
the  no-collision  probability  factor  into  the  model.  Rather,  we  prefer  to 
retain  the  breakdown  surface  model  in  its  present  form  as  a  simple  means  of 
obtaining  an  overestimate  to  the  thermally  backscattered  flux  from  a  centered 
rarefaction  flow.  An  improved  model  can  be  constructed  by  considering  thermal 
backscattering  from  the  entire  flow  field  (tempered  by  the  probability  of 
no-collision),  without  resorting  to  the  physically  untenable  notion  of  an 
abrupt  transition  from  continuum  flow  to  free  molecular  flow. 


5.  CONCLUSIONS 


Some  surprising  similarity  laws  of  the  breakdown  surface  model  were 
observed.  It  has  been  shown  that  they  were  a  direct  result  of  the  self 
similar  structure  of  the  Prandt 1-Meyer  flow  field  to  which  the  model  was 
applied.  Specifically,  it  was  found  (and  shown  plausible)  that  reduced  values 
of  either  the  exhaust  stagnation  density  pQ,  or  the  breakdown  parameter  B, 
caused  higher  backseat tered  flux. 

The  breakdown  surface  model  for  thermally  backscattered  flux  from  a 
centered  rarefaction  fan,  has  been  shown  to  overestimate  the  flux  arriving  at 
the  spacecraft.  It  is  suggested  that  an  improved  model  be  constructed  by 
considering  thermal  backseat  ter ing  from  the  entire  flow  field,  along  with  the 
probability  factor  for  a  side-scattered  molecule  traversing  the  main  flow 
collisionlessly . 

The  molecular  flux  of  corrosive  species  (HF,  DF )  arriving  at  the 
spacecraft  (Figures  6  and  8)  is  no  larger  than  about  107  (sec-1  m-^),  which  is 
negligible  since  it  corresponds  to  about  I0-b  molecular  monolayers  per  year. 
This  conclusion  is  reliable  since  even  this  flux  level  is  an  overestimate. 

The  maximum  thermally  backscattered  flux  of  light  species  (H,  H2,  He)  is 
in  the  range  of  10^u  to  (sec-1  m-11)  (see  Figures  5,  7,  9).  Thus,  we 
conclude  that  while  thermal  backscatter ing  would  contribute  significantly  to 
the  flux  of  light  molecules  arriving  at  the  spacecraft,  it  is  utterly 
negligible  as  far  as  heavy  molecules  are  concerned. 


19 


5.  REFERENCES 


(1)  Chemical  Propulsion  Information  Agency,  "JANNAF  Handbook,  Rocket 

Exhaust  Plume  Technology,  Chapter  6:  Spacecraft  Plume 

Contamination",  CPIA  Publication  263,  June  1983. 

(2)  Chemical  Propulsion  Information  Agency,  "JANNAF  13t'1  Plume  Technology 
Meeting,"  Vol.  I,  Houston,  TX,  April  1982.  CPIA  Publication  357, 
April  1982. 

(3)  G.  A.  Bird,  "Breakdown  of  Continuum  Flow  in  Free  Jets  and  Rocket 
Plumes",  Proc.  12t^1  Symp.  on  Rarefied  Gas  Dynamics.  In  Vol.  74, 
Progress  in  Aeronautics  and  Astronautics,  Sam  S.  Fisher  Editor, 
Published  by  AIAA  1981. 

(4)  H.  B.  Noller,  "Approximate  Calculation  of  Expansion  of  Gas  From 
Nozzles  into  High  Vacuum",  The  Journal  of  Vacuum  Science  and 
Technology,  Vol.  3,  202-207,  1966. 

(5)  F.  Mastrup,  E.  Broadwell,  J.  Miller  and  T.  A.  Jacobs,  "Hydrogen 
Fluoride  Laser  Technology  Study",  Technical  Report  No.  AFWL-TR-72-28 , 
October  1972. 

(6)  G.  A.  Bird,  Molecular  Gas  Dynamics,  Clarendon  Press,  Oxford  1976. 

(7)  G.  A.  Bird,  "Prandtl-Meyer  Flow  of  a  Finite  Knudsen  Number  Gas", 
Proceedings  of  the  Seventh  Australian  Conference  on  Hydrodynamics  and 
Fluid  Mechanics,  1980. 

(8)  R.  J.  Cattolica,  R.  J.  Gallagher,  J.  B.  Anderson  and  L.  Talbot, 

"Aerodynamic  Separation  of  Gases  by  Velocity  Slip  in  Freejet 
Expansions",  AIAA  Journal,  Vol.  17,  P.  344,  1974.  (Also  Reprint 

AIAA-77-709) . 


20 


APPENDIX  A.  The  Computer  Code  RINGBD 


I 

i 

* 

y 


We  present  a  printout  of  the  code  RINGBD  along  with  the  reults  of  the 
nominal  case  (printout  of  actual  run).  This  is  preceded  by  a  brief 
description  of  the  subroutines  and  a  summary  of  major  variables  with  their 
code  and  report  notations. 


A. 1  Description  of  Subroutines 


MAIN  PROGRAM 


INIDAT 


FLUX 


BREAKR 


BREAKM 


-  Computes  flux  integration  by  summation  of  segment 
contributions  (centered).  Printing  of  results. 
Definition  of  all  data  (no  input  file). 
Preparatory  evaluation  of  parameters.  Printing 
of  data. 

-  Evaluates  flux  emitted  from  a  single  point  on 

breakdown  surface  (mean  segment  values)  to  a 
point  on  spacecraft  (XS). 

Computes  point  on  breakdown  surface  for  given 

Mach  number. 

-  Computes  point  on  breakdown  surface  for  mean 

Mach  number  of  a  segment. 


A. 2  Code  Versus  Report  Notation 


XC(I)  -  [AB]  -  Mole  fraction  of  species  AB.  (I«l,2,3,4,5  cor¬ 
responds  to  H,  HF,  H2 »  DF,  He). 

WC(I)  -  Wi  -  Molecular  weight  of  species  i. 

WAV  -  W^  -  Average  molecular  weight 

TO  -  T0  Stagnation  temperature 

RHOO  -  pQ  -  Stagnation  density 

G  -  y  -  Specific  heat  ratio 

EMI  -  M  -  Exit  Mach  number 

LAMDA1-  X^  -  Exit  mean  free  path 
AO  -  Aq  Spacecraft  radius 

R  -  R  -  Distance  from  corner  (X*0,  Y2+Z2  *  A G2). 

DIST  -  L,gs  “  Distance  between  emitting  point  on  breakdown  surface 
and  receiving  point  (XS)  on  spacecraft. 

XS  -  Xg  Point  on  spacecraft  (X-Xg,  Y-A0,  Z-0). 

PSI  -  4i  -  Characteristic  angle 

AMU  -  u  -  Mach  angle 

TETA  -  8  -  Velocity  vector  angle 

PHI  -  4>  -  Rotation  angle  for  flux  integration 

W  -  to  -  Angle  between  x-axis  and  line-of-sight  Lyg 

BETA  -  B  -  Angle  between  Y-axis  and  projection  of  Ljjs  on  (Y,Z) 

plane. 

DM0  -  AM  -  Mach  number  increment  for  flux  integration 

EM  -  M  -  Mach  number 

PBIRD  -  B  -  Breakdown  parameter 


A. 3  Code  Listing  (Run  of  Nominal  Case); 


RINGBD, NOXREF 
IMPLICIT  REAL*8CA-H,0-Z,$) 

COMMON  /GAMA/G,G1,G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G14 
1  G16,G17,G18,G19,G20 

COMMON  /PAR/CO, ENO, EMI, D, TL IM. ETALIM,CLIM, EL 0,Q0, TO, 

1  PBIRD, RBIRD, DM0, DEG, OMEGA, XSV( 51 ) 

COMMON  /NPAR/NETA, NC, NT ,NEMO,NPHI,NXS,NRO, NSPEC 

COMMON  /GEOM/APF ,PAI,PAI2,SW,CW,BETA,SBETA,CBETA,PSI1,SPSI1, 

1  CPSI1,PSIF,SPSIF,CPSIF,AK,SK,CK,A0,RF,XF,YF,ZF, 

2  PHI,SPHI,CPHI, RMIN, RMAX, XS , DIST, 

3  AMU1 ,ZETA1,XN,YN,ZN,PSIM,SPSIM,CPSIM,R0 
COMMON  /EPSIL/EPSQ, EPSETA, EPST, EPSC, EPSEM 

COMMON  /EXTREM/TEXT, ETA EXT , CEXT, REXT , PSI EXT, EMEXT, BEXT, OEXT 
COMMON  /SPEC/ WAV , XC(5),WC(5),WCR(5), XNAMEC 5) , QCC 5 ) , FLUXC( 5 ) 
DIMENSION  DSUMC5) 

PRINT  101 

FORMAT C'l'/lX, • RINGBD  -  FLUX  INTEGRATION  FROM  BREAKDOWN', 

1  IX, 'SURFACE'//) 


CALL  INIDAT 
PRINT  110,XNAME 

FORMAT (///IX, '  NX','  NEM','  XS  ',' 

1  5( 4X, A6 , IX, '/  LOG', IX)) 

DO  200  NX=1 , NXS 
EM=EM1 

CALL  BREAKR(EM, RF) 

IF(NRO.GT.O)  RF=RO 
XF=RF*CPSI1 
YF=RF*SPSI1+A0 
XS=XSV(NX) 

QMAX=0 . 

DO  45  N=l, NSPEC 
FLUXC(N)=0 . 

CONTINUE 

DO  1  NEM=1 , NEMO 

RN=RF 

XN=XF 

YN=YF 

IF(NEM.GT.NRO)  GO  TO  41 

RF=RO+DFLOAT(NEM)X(RMIN-RO)/DFLOAT(NRO) 

XF=RF*CPSI1 

YF=RF*SPSI1+A0 

RMEAN=( RN+RF)/2 . DO 

EMMEAN=EM1 

PSIM=PSI 1 

SPSIM=SPSI1 

CPSIM=CPSI 1 

GO  TO  42 

CONTINUE 

EM=EM+DM0 

CALL  BREAKKf CM, RF) 

EMMEAN=EM-DM0/2 . DO 
CALL  BREAKM( EMMEAN,RMEAN) 

CONTINUE 

AL0NG=DSQRT((XF-XN)X*2+(YF-YN)XX2) 
SALFA=(YF-YN)/ALONG 
CAL FA=(XF-XN)/A LONG 
PHIMAX=DARCOS( AO/ ( AO+RMEANXSPSIM) ) 

DPHI =PHIMAX/NPHI 


PHIMAX', 


RIN00010 
RIN00020 
,  G15 , RINOO 0 30 
RIN00040 
RIM 00050 
RIN00060 
RIN00070 
RIN00080 
RIN00090 
RIN00100 
RIN00110 
RIN00120 
RIN00130 
RIN00140 
RIN00150 
RIN:  160 
RIN00170 
RIN00180 
RIN00190 
RIN00200 
RIN0021 0 
RIN00220 
',  RIN00230 

RIN00240 
RIN00250 
RIN00260 
RIN0027  0 
RIN00280 
RIN00290 
RIN00300 
RIN0031 0 
RIN00320 
RIN00330 
RIN00340 
RIN00350 
RIN00360 
RIN0037  0 
RIN00380 
RIN00390 
RIN00400 
RIN0041 0 
RIN00420 
RIN00430 
RIN00440 
RIN00450 
RIN0046  0 
RIN0047  0 
RIN00480 
RIN00490 
RIN00500 
RIN0051  0 
RIN00520 
RIN00530 
RIN00540 
R1N00550 
RIN00560 
RIN0057  0 
RIN00580 
RIN 00590 
RIN00600 
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DO  44  N=1 , NSPEC 
DSUM( N ) =0 . 

CONTINUE 

DO  2  NP  =  1 , NPHI 

PHI  =  <  DFLOAT (NP)-0.5D0)XDPHI 

CALL  FLUX( EMMEAN, RMEAN) 

CR0SS1=SW*CBETA 

CR0SS2=(SALFA)»(-CW)+(-CALFA*CPHI)X(-SWXCBETA)+ 

1  (-CALFAxSPHI)x(-SWxSBETA) 

GOREM=CROSS1XCROSS2*DPHI*(AO+RMEANXSPSIM)XALONG/DIST*X2 
DO  24  N=l, NSPEC 
DSUM(N)=DSUM(N)+QC(N)XGOREM 
CONTINUE 

IF(QMAX.GE.QEXT)  GO  TO  25 

QMAX=QEXT 

CONTINUE 

PRINT  22, NEM,NP, EMMEAN, RMEAN, PHIMAXXDEG,DARCOSCCW)XDEG, 

1  BETAXDEG,PHIXDEG,CR0SS1,CR0SS2, ALONG, DIST, GOREM, 

2  SPSIM, AO, ( QCCN), FLUXCC N),N=1, NSPEC) 

F0RMAT(/1X, ' NEM, NP , EMMEAN , RMEAN, PHIMAX= ',2I3,3D13.4/ 

1  1X,'W,BETA,PHI=',3D15.5/ 

2  IX, ' CR0SS1 , CR0SS2, ALONG, DIST, GOREM= ',5015.5/ 

3  IX, 'SPSIM, A0=' ,2015.5/ 

4  IX, ' QC, FLUXC=  *,5(1X,D10.3,1X,D10.3)) 

CONTINUE 

DO  26  N=l, NSPEC 
FLUXC(N)=FLUXC(N)+DSUM(N) 

CONTINUE 

IFCNEM.LE.NR0+2)  GO  TO  1 
DO  27  N=l, NSPEC 

IF(CDSUM(N)/FLUXC(N)) .GT.EPSEM)  GO  TO  28 

CONTINUE 

GO  TO  10 

CONTINUE 

CONTINUE 

CONTINUE 

DO  31  N=l, NSPEC 

FLUXC(N)=2 . DOxXCCN)xFLUXC(N) 

CONTINUE 

PRINT  11 , NX, NEM, XS, PHIMAXxDEG.QMAX, 

1  ( FLUXCC  N ) , DABS ( DL0G1 0 ( FLUXCC  N))),N  =  1, NSPEC) 

FORMAT (/1X,2I4,F9.4,F7 .2, Dll .3,5(1X, DIO. 3, ’/»,F5.2)) 

CONTINUE 
PRINT  102 

F0RMAT(///1X, 'END  RINGBD  RUN',///) 

STOP 

END 

SUBROUTINE  INIDAT 
IMPLICIT  REALX8(A-H,0-Z,$) 

REALX8  LAMDAO, LAMDA1 

COMMON  /GAMA/G, G1 ,G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G14, 
1  G16 , G17 ,G18,G19,G20 

COMMON  /PAR/CO, ENO, EMI ,D,TLIM, ETALIM.CLIM, ELO,QO,TO, 

1  PBIRD, RBIRD, DMO, DEG, OMEGA, XSV< 51 ) 

COMMON  /NPAR/NETA, NC, NT , NEMO, NPHI , NXS , NRO , NSPEC 

COMMON  /GEOM/APF, PAI,PAI2,SW,CW,BETA,SBETA,CBETA,PSI1,SPSI1, 

1  CPSI1,PSIF,SPSIF,CPSIF,AK,SK,CK,A0,RF,XF,YF,ZF, 

2  PHI, SPHI,CPHI,RMIN,RMAX,XS, DIST, 

3  AMU1 ,ZETA1,XN,YN,ZN,PSIM, SPSIM, CPSIM,RO 
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95 

COMMON  /EPSIL/EPSQ, EPSETA, EPST, EPSC, EPSEM 

RIN01200 

96 

COMMON  /EXT REM/ TEXT , ETA EXT, C EXT, R EXT , PS I EXT , 

EMEXT , BEXT , QEXT 

RIN01210 

97 

COMMON  /SPEC/ WAV , XC(5),WC(5),WCR(5),XNAME(5) 

, QC( 5 ) , FLUXCC  5 ) 

RIN01220 

98 

DATA  XC/.091D0, .091D0, .104D0, .135D0, .579D0/ 

RIN01230 

99 

DATA  WC/1 .00DO,20.0DO,2.0OD0,21 .0D0,4. OODO/ 

RIN01240 

100 

DATA  XNAME/'  H  HF  H2  DF 

HE  '/ 

RIN01250 

101 

PAI=4 . KDATANC . ID  1) 

RIN01260 

102 

AR=8 . 3143D3 

RIN0127  0 

103 

AV=6 . 022D  26 

RIN01280 

C 

OMEGA=0 . 5  IS  FOR  HARD  SPHERE  COLLISIONS, 

RIN01290 

C 

AN 

AVERAGE  RECOMMENDED  VALUE  IS  ABOUT  0MEGA=0.75 

RIN01300 

104 

OMEGA=0 . 5D0 

RIN01310 

105 

NSPEC=5 

RIN01320 

106 

WAV=0. 

RIN01330 

107 

DO  51  N=1 , NSPEC 

RIN01340 

108 

WAV=WAV+XC(N)XWC(N) 

RIN01350 

109 

51 

CONTINUE 

RIN01360 

110 

DO  52  N=l, NSPEC 

RIN01370 

111 

WCR(N)=DSQRTCWC(N)/WAV) 

R I  NO  1380 

112 

52 

CONTINUE 

RIN01390 

113 

A0=2 . 5D0 

RIN01400 

114 

EMI =4 . 0D0 

RIN01410 

115 

RHO0=0 . 0075D0 

RIN01420 

116 

T0=2 . 300D3 

RIN01430 

117 

D=2 . 50D-10 

RIN01440 

118 

G=1 . 54D0 

RIN01450 

119 

ENO=RHOO*AV/WAV 

RIN01460 

120 

C0=DSQRT(G*AR*T0/WAV) 

RIN0147  0 

121 

PBIRD=0.05D0»2.D0 

RIN01480 

C 

R0 

IS  THE  RADIUS  FOR  BEGINNING  THE  INTEGRATION 

ALONG  THE  M=M1 

RIN01490 

C 

CHARACTERI STIC ( THE  AUGMENTED  BREAKDOWN  SURFACE) 

RIN01500 

C 

NR0  IS  THE  NUMBER  OF  INTEGRATION  INTERVALS  ON  THIS  SEGMENT. 

RIN01510 

C 

FOR  NO  INTEGRATION  ALONG  M=M1  CHARACTERISTIC,  SET  NR0=0. 

RIN01 520 

122 

R0  =  0 . 

RIN01530 

123 

NR0=10 

RIN01540 

124 

DM0=0 . IDO 

RIN01550 

125 

NEM0=20 . D0/DM0+NR0 

RIN01560 

c 

TO 

GET  FLUX  DUE  TO  AUGMENTED  BREAKDOWN  SURFACE 

SOLELY,  ACTIVATE! 

RIN0157  0 

c 

NEM0=NR0 

RIN01580 

126 

NPHI = 1 0 

RIN01590 

127 

NXS=13 

RIN01600 

128 

XSI=1 . D-2 

R1N01610 

129 

XSF=1 . D1 

RIN01620 

130 

XSV(1)=XSI 

RIN01630 

131 

IF(NXS.EQ.l)  GO  TO  111 

RIN 01640 

132 

DXL=(DLOG(XSF)-DLOG(XSI) )/( NXS-1 . DO ) 

RIN01650 

133 

XLI=DLOG(XSI) 

RIN01660 

134 

DO  11  NX=2 , NXS 

RIN01670 

135 

XSVC  NX ) =DEXP(XL I  +  CNX-1 . DO ) XDXL ) 

RIN01680 

136 

11 

CONTINUE 

RIN01690 

137 

111 

CONTINUE 

RIN01700 

138 

EPSEM=1 . D-5 

RIN01710 

139 

DEG=180 . DO/PAI 

RIN01720 

140 

PAI2=PAI/2 . DO 

RIN 01730 

141 

GAMMA=G 

RIN01740 

142 

G1  =  (G-1 . DO )/2 . DO 

RIN017  50 

143 

G2=(G+1 .D0)/(2.D0X(G-1 .DO)) 

RIN01760 

144 

G3=G/2 . DO 

RIN01770 

145 

G4=(G+1 . DO )/(G-l . DO ) 

RIN01780 

146 

G5=DSQRT ( (G+l . D0)/(G-1 . DO) ) 

RIN01790 

a 
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147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 


176 


177 


178 

179 

180 
181 


G6=l . D0/(G-1 . DO ) 

G7=2 . D0/(G+1 . DO ) 

G8=(5.D-1X(G+1.D0)XX2/(G-1.D0))XX(1 . DQ/CG+l .DO))  x 
1  (CG+1.DO)/(G-1.DO))XX(CG-1.DO)/(G+1.DO)) 

G9=CG+3.D0)/(2.D0XCG-1.D0)) 

G10=(7 .DO-3.DOXG)/(2.DOX(G-1 .DO)) 

G11=DSQRT(G/PAI)/(2.D0X(G+1 .DO)) 

G12=DSQRT(G/2.D0) 

G13  =  l  .D0/DSQRK2.D0XGXPAIXX3) 

LAMDA0=1 .D0/(DSQRT(2.D0)XPAIXDXX2XEN0) 

LAMDA1=LAMDA0XU . DO+G1XEM1XX2)XX(G6-OMEGA+Q . 5D0) 

RBIRD=G11/(DXX2XEN0XPBIRD) 

ZETA1=G5*DATAN(DSQRTCEM1XX2-1.D0)/G5) 

AMU1=DARSIN( 1 . D0/EM1 ) 

P5I1=PAI2+AMU1 

SPSI1=DSIN(PSI1) 

CPSI1=DC0S(PSI1) 

PSIF=PAI2+AMU1+ZETA1-G5XPAI2 

SPSIF=DSIN(PSIF) 

CPSI F=DCOS( PSIF ) 

CALL  BREAKR(EMl.RMIN) 

RSMIN=RMINXDSQRT( C2.D0/(G-1 .D0)+EM1XX2)/(EM1XX2-1 .DO)) 

C 

PRINT  201 # NSPEC, XNAME 

201  FORMAT(/lX, 'SPECIES  DATA  NSPEC=',I3/ 

1  IX, 'SPECIES  NAMES  1 , 11 (2X, A6 , 2X) ) 

PRINT  202, XC 

202  FORMAT (  IX, 'MOL E  FRACTION  XC= ', 1 1 C F8 . 4, 2X) ) 

PRINT  203, WC 

203  FORMAT (  IX, 'MOL.  WEIGHT  WC= * , 1 1 C F8 . 4, 2X) ) 

PRINT  21,AR,AV,WAV,G,RHOO,TO,ENO,CO,D 

21  FORMATC/1X, 'THERMODYNAMIC  DATA'/ 

1  IX, ’ AR , AV, WAV , GAMMA = ' , 2X, 2D14 . 5, 2F9 . 3/ 

2  IX, ’RHO0,T0,EN0,C0,D=*,D12.4,F8.0,D13.5,2D12.4) 

PRINT  22, EMI , PSI1*DEG, PSIFxDEG,PBIRD, 

1  AO , RMIN, RSMIN, RO , 

2  LAMDAO , LAMDA1 , RMIN/LAMDA1 , RSMIN/L AMDA1 

22  FORMATC/1X, 'FLOW  AND  GEOMETRY  DATA'/ 

1  IX, 'EM1,PSI1,PSIF,PBIRD=' ,4F9.3/ 

2  IX,  'A0,RMIN,RSMIN,R0=',F9.3,3D14.5/' 

3  IX, 'LAMDAO , LAMDA1 , RMIN/LAMDA1 , RSMIN/LAMDA1= ' , 4D13 . 4) 
PRINT  23, DM0 , NPHI , NRO , EPS EM 

23  F0RMAT(/1X, 'INTEGRATION  DATA'/ 

1  IX, ' DM0 , NPHI , NRO , EPSEM= ',F9.3,2I5,4X,D12.3) 

RETURN 

END 
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SUBROUTINE  FLUX(EM,R) 

IMPLICIT  REALX8(A-H,0-Z,$) 

COMMON  /GAMA/G , G1,G2,G3,G4,G5,G6,G7,G8,G9,G10,G11,G12,G13,G14 
1  G16,G17,G18,G19,G20 

COMMON  /PAR/CO, ENO, EMI ,D,TL IM,ET ALIM, CL IM, EL0,Q0,T0, 

1  PBIRD,RBIRD, DM0 , DEG, OMEGA , XSV ( 51 ) 

COMMON  /NP AR/ NET A , NC, NT , NEMO , NPHI , NXS, NRO , NSPEC 

COMMON  /GEOM/APF, PAI,PAI2,SW,CW, BET A , SBETA , CBETA , PSI 1 , SPSI 1 , 

1  CPSI1 , PSI F, SPSI F.CPSIF, AK,SK,CK,AO,RF,XF,YF,ZF, 

2  PHI, SPHI.CPHI, RMIN, RMAX , XS , DI ST , 

3  AMU1 ,ZETA1,XN,YN,ZN,PSIM,SPSIM,CPSIM,R0 
COMMON  /EXTREM/TEXT , ETAEXT , CEXT , REXT , PSI EXT , EMEXT , BEXT , QEXT 
COMMON  /SPEC/ WAV , XC( 5 ) , WC( 5 ) , WCRC  5 ) , XNAME ( 5 ) , QC( 5 ) , FLUXC( 5 ) 


RIN0226  0 
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191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
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224 

225 

226 

227 

228 


SPHI=DSIN(PHI ) 

CPHI=DC0SCPHI) 

XMEAN=RXCPSIM 
YMEAN=(AO+R*SPSIM)*CPHI 
ZMEAN=( AO+RXSPSIM ) xSPHI 
TBETA=ZMEAN/( YMEAN-AO ) 

BETA=PAI2-DATAN( 1 . DO/ T BETA) 

SBETA=DSIN(BETA) 

CBETA=DCOS(BETA) 

DIST=DSQRT( ( XS-XMEAN) *X2+( YMEAN-AO )XX2+ZMEANX*2) 

CW=-(XS-XMEAN)/DIST 

SW=DSQRT ( 1 . D0-CWXX2) 

GM=(1 ,D0+G1XEMX*2)X*(-G2) 

AMU=DARSIN( 1 . DO/EM) 

TETA=PSIM-AMU 

STETA=DSIN(TETA) 

CTETA=DCOS(TETA) 

CKAPA=(CTETA)X(-CW)+(STETAXCPHI)X(-SWXCBETA)+ 

1  (STETAXSPHI )X(-SWXSBETA) 

SKAPA=DSQRT(1 .D0-CKAPAXX2) 

QEXT=0 . 

DO  1  N  =  1 , NSPEC 
EMT=EMXCKAPAXG12XWCR(N) 

IF(DABSCEMT) .GT. 13. DO)  EMT=EMTX( 1 3 . DO/DABSC EMT) ) 

P0W=G3XEMXX2XWCR(N)XX2 

POWT  =  P0WXSKAPAX*2 

IFCPOW.GT.l .OD2)POW=1 .0D2 

I F ( POWT . GT . 1 . 0D2 )  POWT=1.0D2 

EXP1 =DEXP( -POW) 

EXP2=DEXP( -POWT ) 

ERFC1 =DERFC( -EMT ) 

IF(ERFCl.LT.l.D-43)  ERFCl=l.D-43 
C  IF(XS.LT.5.D0)  GO  TO  234 

C  IF( EXP1 . GT . 1 . D-20 . AND . ERFC1 . GT . 1 . D-20 )  GO  TO  234 
C  PRINT  235, N,POW, POWT, EMT, EXP1 , EXP2, ERFC1 

C235  FORMAT ( I4,8D11.3) 

C234  CONTINUE 

EVER1=EN0X(C0/WCR(N))XG13XGMX(1 ,D0+EMTxx2)*EXPl 
EVER2=( ENOx( CO/WCRC  N ) ) xo . 5D0/PAI ) x( EMxWCRI N) ) xCKAPAxGMx 
1  (1 ,5D0+EMTxx2)XEXP2XERFCl 

QC(N)=EVER1+EVER2 
IF( QEXT . GE . QC( N ) )  GO  TO  1 
OEXT=OC(N) 

1  CONTINUE 

RETURN 
END 


RIN02390 
RIN 02400 
RIN02410 
RIN02420 
RIN02430 
RIN02440 
RIN 02450 
RIN02460 
RIN0247  0 
RIN02480 
RIN02490 
RIN 02500 
RIN0251 0 
RIN02520 
RIN02530 
RIN02540 
RIN02550 
RIN02560 
RIN0257  0 
RIN02580 
RIN02590 
RIN02600 
RIN02610 
RIN02620 
RIN02630 
RIN02640 
RIN02650 
RIN0266  0 
RIN0267  0 
RIN02680 
RIN02690 
RIN027  00 
RIN027 1 0 
RIN02720 
RIN027  30 
RIN02740 
RIN02750 
RIN027  6  0 
RIN0277  0 
RIN02780 
RIN02790 
RIN02800 
RIN02810 
RIN02820 
RIN02830 
RIN02840 


229  SUBROUTINE  BREAKR(EM,R)  RIN02850 

230  IMPLICIT  REALX8(A-H,0-Z,$)  RIN02860 

231  COMMON  /GEQM/ APF , PAI , PAI2, SW, CW, BETA, SBETA , CBETA , PSI 1 , SPSI 1 ,  RIN02870 

1  CPSI1,PSIF,SPSIF,CPSIF,AK,SK,CK,A0,RF,XF,YF,ZF,  RIN02880 

2  PHI , SPH3,CPHI , RMIN, RMAX, XS , DIST ,  RIN02890 


3  AMU1 ,ZETA1,XN,YN,ZN,PSIM,SPSIM,CPSIM,R0  RIN02900 

232  COMMON  /PAR/CO , ENO , EMI , D , TL IM, ETAL IM, CL IM, EL  0 , QO , TO ,  RIN02910 

1  PBIRD,RBIRD, DMO, DEG, OMEGA, XSVL51)  RIN02920 

233  COMMON  /GAMA/G, Cl , G2 , G3 , G4 , G5, G6 , G7 , G8, G9 , G10, Gil ,G12, G13, G14 , G15 , RIN02930 

1  G16,G17,G18,G19,G20  RIN02940 

234  R=RBIRDXDSQRT ( EMXX2-1 .D0)X(1 . D0+G1XEM*X2 )*x(G6-0MEGA+0 . 5D0 )  RIN02950 

235  ZETA=G5XDATAN(DSQRT(EMXX2-1.D0)/G5)  RIN02960 

236  PSI=PAI2+AMU1+ZETA1-ZETA  RIN02970 


237  XF=R*DCOS( PSI ' 

238  YF=R*DSIN(PSI)+AO 

239  ZF=0. 

240  1  CONTINUE 

241  RETURN 
C 

242  ENTRY  BREAKM(EM,R) 

243  C  R=RFIRD*DSQRT(EM**2-1 . DO )»( 1 . D0+G1*EM**2)**(G6-OMEGA+0 . 5D0) 

244  ZETA=G5*DAT  AN(DSQRT(EM**2-1.D0)/G5) 

245  PSIM=PAI2+AMU1 +ZETA1-ZETA 

246  SPSIM=DSIN(PSIM) 

247  CP$IM=DCOS(PSIM) 

248  RETURN 

249  END 


RIN02980 
RIN02990 
RIN03000 
RIN 03010 
RIN03020 
RIN 03030 
RIN03040 
RIN 03050 
RIN03060 
RIN 03070 
RIN03080 
RIN 03090 
RIN 03100 
RIN03110 
RIN03120 


SENTRY 


RIN03130 


-  FLUX  INTEGRATION  FROM  BREAKDOWN  SURFACE 


SPECIES  DATA  NSPEC«  5 

SPECIES  NAMES  H  HF  N 2  DF  MF 

MOLE  FRACTION  XC«  0.0910  0.0910  0.1090  0.1350  0.5790 

MOL.  HEIGHT  NC*  1.0000  20.0000  2.0000  21.0000  9.0000 

THERMODYNAMIC  DATA 

AR.AV.WAV, GAMMA*  0.83143D04  0.60220027  7.270  1.540 

RHOO.TO.ENQ.CO.D*  0.7500D-02  2300.  0.62125D  24  0.2013D  04  0.2500D-09 

FLOW  AND  GEOMETRY  DATA 

EM1.PSU.PSIF.P81RD*  4.000  104.478  41.044  0.100 

AO.RMIN.RSMIN.RO*  2.500  0  50374D-02  0.34812D-02  0.00000D  00 

wAMOAO.LAMDAl.RMlN/LAMDAl.RSMIN/LAMDAl>  0.5797D-OS  0. 12810-03  0.2372D  02  0.2718D  02 

INTEGRATION  DATA 

DM0 , NPHI . NRO , EPSEH*  0.100  10  10  P.100D-04 


NX 

NEM 

XS 

PHIMAX 

QMAX 

H 

✓  LOO 

HF 

/ 

LOG 

M2 

/  LOO 

DF 

/ 

LOO 

HE 

/  LOO 

1 

47 

0.0100 

11.19 

0.120D 

25 

0 . 8240 

22/21.92 

0.2690 

07/ 

6.43 

0.7770 

21/20.89 

0.691D 

06/ 

5.84 

0.6530 

20/19.81 

2 

50 

0.0178 

12.15 

0.1210 

25 

0.474D 

22/21 .68 

0.1620 

07/ 

6.21 

0.424D 

21/20. *3 

0.4180 

06/ 

5.62 

0.3490 

20/19.54 

3 

54 

0.0316 

13.50 

0.1220 

25 

0.2370 

22/21.38 

0.7760 

06/ 

5.69 

0.1960 

21/20 . 29 

0.202D 

06/ 

5.30 

0.1540 

20/19.19 

4 

59 

0.0562 

15.27 

0.123D 

25 

0. 106D 

22/21.02 

0 .2950 

06/ 

5.47 

0.7800 

20/19.89 

0.7700 

05/ 

4.89 

0.568D 

19/18.75 

5 

64 

0.1000 

17.16 

0.1240 

25 

0.421D 

21/20.62 

0.892D 

0  5/ 

4.95 

0.268D 

20/19.43 

0.233D 

05/ 

4.37 

0.175D 

19/18.24 

4 

70 

0.1778 

19.55 

0. 1240 

25 

0.149D 

21/20.17 

0.219D 

0  5/ 

4 . 34 

0.7910 

19/18.90 

0.5750 

04/ 

3.76 

0.4510 

18/17.65 

7 

77 

0.3162 

22.52 

0. 124D 

25 

0 .471  D 

20/19.67 

0.464D 

04/ 

3.67 

0.2030 

19/18.31 

0. 1220 

04/ 

3.09 

0 .1000 

18/17.00 

8 

84 

0.5623 

25.65 

0.124D 

25 

0.132D 

20/19.12 

0.9010 

03/ 

2.95 

0.4610 

18/17.66 

0.2370 

03/ 

2.37 

0.201D 

17/16.30 

9 

92 

1.0000 

29.39 

0. 1 24D 

25 

0.335D 

19/28.52 

0.167D 

03/ 

2.22 

0.959 D 

17/16 . 98 

0.4400 

02/ 

1  .64 

0. 381D 

16/15.58 

10 

100 

1.7783 

33.26 

0.124D 

25 

0.771D 

18/17.89 

0 . 305D 

02/ 

1.48 

0. 167D 

17/16.27 

0  6010 

01/ 

0.90 

0 .6990 

15/14  84 

11 

108 

3.1623 

37.22 

0  .  1  24  D 

25 

0.1650 

18/17.22 

0 . 548  D 

01/ 

0.74 

0.3530 

16/15.55 

0. 1440 

01/ 

0.16 

0. 126D 

15/14.10 

12 

116 

5.6234 

41.19 

0. 124D 

25 

0.3320 

17/16.52 

0  9610 

00/ 

0.01 

0.6480 

15/14.81 

0.2580 

00/ 

0 . 59 

0.227D 

14/13  36 

13 

124 

10.0000 

45.14 

0.124D 

25 

0.641D 

16/15.81 

0  1750 

00/ 

0.76 

0.1170 

15/14  07 

0.4610- 

01/ 

l  .  34 

0.405D 

13/12  61 

END  RlNGBO  RUN 
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Thermal  Backseat 
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FLUX  OF  SPECIES  HF 


[hf]  =  .091 


Figure  6.  Flux  of  Species  HF  at  Various  Stagnation  Densities 


I 


34 


.s.nV..v-.V^ 


FLUX  (molecules  sec* 


1022! 


FLUX  OF  SPECIES  H2  [h2]  =  .104 

M-j  =4  (f=  1.54 

T0  =  2300  (K)  WA  =  7.27  B  =  .05 


Figure  7.  Flux  of  Species  H2  at  Various  Stagnation  Densities 
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Figure  9.  Flux  of  Species  He  at  Various  Stagnation  Densities 
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